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Abstract. We consider bipartite quantum systems that are described completely by a state vector ^(t)} 
and the fully deterministic Schrodinger equation. Under weak constraints and without any artificially in- 
troduced decoherence or irreversibility, the smaller of the two subsystems shows thermodynamic behaviour 
like relaxation into an equilibrium, maximization of entropy and the emergence of the Boltzmann energy 
distribution. This generic behaviour results from entanglement. 

PACS. 05.70.Ln Nonequilibrium and irreversible thermodynamics - 05. 30.-d Quantum statistical mechan- 
ics 



1 Introduction 

Various attempts have been made for quantum systems to 
account for thermodynamic behaviour pQ [2] 0, in par- 
ticular for the relaxation into an equilibrium state. Anal- 
ogous to the Hamilton equations in classical mechanics, 
the microscopic equation determining quantum mechani- 
cal evolution, the Schrodinger equation, is time-reversible, 
i.e. invariant under time reversal. Therefore in quantum 
mechanics there exist the same problems as in classical 
mechanics 0] when one wants to deduce irreversible be- 
haviour of thermodynamic systems. 

One way to 'introduce' irreversibility in quantum me- 
chanics is the treatment of open systems by means of 
master equations [5] In this paper we consider an- 
other very general approach to this problem that does 
not need those specific assumption about an environment 
that is not part of the system considered. Instead, sys- 
tem and environment arc treated as a whole via pure 
Schrodinger dynamics; even for small systems (i.e. small 
accessible Hilbert-spaces) thermodynamic behaviour and 
quasi-irreversibility is shown to emerge for a large class of 
systems. 

We have performed numerical simulations on relatively 
small quantum systems (up to total state vector dimen- 
sions of N tot w 1000). They are partitioned into a smaller 
"part of interest", the 'gas' g, and a larger part, the envi- 
ronment or 'container' c, N tot = N 9 N C . The interaction 
between those two subsystems has to be small in com- 
parison to the local energies (i.e. (H mt ) <c {{H 9 ), (H c )} 
for typical states). This constraint as well as the parti- 
tioning is typical for all thermodynamic systems. Even if 
there is no energy exchanged between the two subsystems 



(microcanonical coupling), an environment ('container') is 
eventually needed to maintain the extensive parameters of 
the system considered like e.g. volume etc. As opposed to 
the open system approaches based on a bath, our envi- 
ronment c does have neither a temperature nor a specific 
energy distribution. Its main function is to allow for en- 
tanglement between the two subsystems c and g. 

The numerical simulations shown in this paper are in- 
tended to test wether this entanglement may, indeed, be 
considered the origin of thermodynamics and irreversibil- 

ity mmm 



2 The model 

The Hamiltonian of the full system will consist of the sum 
of the two local Hamiltonians of the gas and the container, 
respectively, plus the interaction Hamiltonian: 



H = H<® 



H mt = H 9 ® l c + I 9 <g> H c 



(1) 



Here, H 9 includes already the confining effect of the con- 
tainer by means of an effective potential; H mt describes 
dynamical corrections (see [S]). The present analysis will 
focus on 'typical' Hamiltonians rather than on concrete 
physical models based on specifically interacting particles 
(as e.g. in JO] |H|k Our system will thus be defined in 
terms of abstract energy spectra and statistical interac- 
tions. 

We consider the eigenbasis of H(°\ the product-states 
\i) 9 ® |m, s) c with i = 0, 1, 2, . . . , N 9 - 1 denoting the 
non-degenerate energy eigenstates of g, m = 0,1,2,... 
the energy eigenstates of c with corresponding degeneracy 
s = 1, 2, . . . , n c m . Any state vector can thus be written as 
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The Hamiltonian models will be assumed to be defined 
with respect to this basis. is thus diagonal. For most 
of the simulations, the simplest model for the gas is used: 
a two-level system with level spacing AE 9 defining the 
pertinent energy scale. 

We assume that the product state basis is in no way 
a preferred basis for the interaction H mt . Its matrix is 
therefore treated as a hermitian random matrix, spec- 
ified by some distribution function u^if?™*). If nothing- 
was known, this distribution function should not depend 
on the representation chosen, i.e. must be invariant un- 
der any unitary transformation. This is guaranteed by 
taking oc exp[— jTr{(H mt ) 2 }] implying indepen- 

dent zero-mean Gaussian distributions with variance <7o = 
^/2AE 9 a for the diagonal matrix elements and variance 
a = AE 9 a for the real and imaginary parts of the off- 
diagonal elements We induce weak coupling by re- 
quiring a <C 1, i.e. the energy scale of H mt is small com- 
pared with that of H^°\ Microcanonical constraints are 
built in by setting all blocks within H mt equal to zero, 
which would connect different energy states within g. 
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Fig. 1. Pertinent relevant energy spectrum for the case of 
microcanonical coupling. 

the subsystem as well as for the entanglement between the 
two subsystems: 



S 9 {p 9 ) = -k B Tr{p 9 lnp 9 } 



(5) 



(with the Boltzmann-constant ks)- S 9 is zero initially 
(product state) but not conserved. 

In our simulations we set h = 1 and AE 9 = 1 which 
gives time the unit of . 



3 Initial state and simulation 

To observe relaxation one has to choose rather special ini- 
tial states, since by far most of the states in the acces- 
sible region of total Hilbert space are equilibrium states 
with respect to the smaller subsystem g jS] HJ. In ^U] e.g. 
randomly chosen states are taken as initial states and no 
relaxation can be observed. Here we choose product states 
of the two subsystem states \ifj 9 ) and \ip c ): 



\&{t = o)) = \r) ® iv> c >- 



(3) 



With this choice we also avoid accidentally taking an eigen- 
state of the total Hamiltonian H (including H mt ) which 
would lead to a plain oscillation but not to a decay. We as- 
sume to know the corresponding energy distribution p tot (E 
which is a constant of motion (closed total system). From 
ptot we can ca i cu i a te the mean energy U — E tot . 

The state vector \&{t)) of the total system will be cal- 
culated numerically using a diagonalisation of the total 
Hamiltonian H. Through H mt the two subsystems g and 
c entangle and can no longer be described by local state 
vectors \ip 9 ) and \ip°). Local properties of the subsytems 
are described using the reduced density matrix of e.g. the 
gas: 



EE< 



4 Results 

4.1 Microcanonical coupling 

The microcanonical coupling in thermodynamics excludes 
energy exchange between the gas and the container, i.e. 
the energy distribution of the subsystems, p 9 and p c , are 
separate constants of motion. Therefore, if we restrict our- 
selves to sharp energy initial states in the container, only 
one energylevel with degeneracy N c has to be considered 
in this subsystem. Fig. ^ shows the energy spectrum for 
the case of microcanonical coupling. Eq and Ef are the 
energies of the two levels in the gas, E c of the container 
level. For a system as shown in fig. ^ with N c = 50 
(i = 0,1; m = Q;s = 1,2,...7V C ), fig. H shows the be- 
haviour of a/|Pqi| 2 (the off-diagonal element of the re- 
duced density matrix of the 2-level-gas) and S 9 (its von- 
*°*)^eumann entropy) . In a system with microcanonical cou- 
pling, the diagonal elements of the reduced density ma- 
trix (here p 9 m = 0.15; p\ x = 0.85) are constants of mo- 
tion, but the off-diagonal clement decreases from the value 
of the product state V0.15 • 0.85 = 0.3571 towards zero, 
the value for a completely mixed state. The average en- 
tropy in equilibrium of this simulation (after ten units of 
time) is (in units of ks) S 9 w 0.4157, the maximal en- 
tropy Sljbaz. ~ 0.4226. The largest possible value for this 
setting would be that of a completely mixed state with 
Pq! = and can be calculated using the eigenvalues of 



i; to, s 



\p\jv> 



(4) P J 



S 



rnax^th 



0.4227 which agrees with the results from 



where p=\W)(\P\is the density matrix of the total system. 
The diagonal elements of the reduced density matrix of 
the gas, pfj, are the probabilities of finding the gas in the 
corresponding local state \i) 9 . The offdiagonal elements of 
p 9 are a measure for local coherence. 

From eq. fllf . the reduced von-Neumann entropy S 9 
of the gas can be calculated, a measure for the purity of 



the simulation. In all the simulations done, S 9 is always 
smaller than S^L th and the difference between those two 
values decreases with increasing system size N c . 



4.2 Canonical coupling 

We still restrict ourselves to a sharp initial energy in c, Ef. 
Conservation of the total energy distribution requires the 
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Fig. 2. '2-level-gas' under microcanonical constraints: Von- 
Neumann entropy S 9 (in units of fes) and absolute value of 
the off-diagonal element of the reduced density matrix as a 



function of time in units of 
N c = 50. 
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Fig. 3. Pertinent relevant energy spectrum for the case of 
canonical coupling and sharp initial energies E\ in the con- 
tainer, rij axe the respective degeneracies. 



consideration of only two more energies: Eq—E\ — AE 9 
and E% = Ef + AE 9 . Fig. [3] shows this energy spectrum. 
Each of the three energylevels of the container m = 0, 1, 2 
is degenerate and contains n c m levels. 

The spectrum of the container in between those three 
levels does not play any role. Important, however, are the 
corresponding degeneracies. Here the usual assumptions 
about large systems are made ^3], i- e - the degeneracy of 
the energy levels is taken to grow exponentially with the 
energy: n c (E c ) oc e yE , which we call the thermodynamic 
scheme of degeneracy. 

Fig. 01 shows the probability p^ Q of finding the sub- 
system g in the ground state |0) ff plotted over time for 
three different simulations that all had the same energy 
spectrum: Eq = 0, Ef = 1 (both non-degenerate); Eq = 
0;n c o (E c Q ) = 50; Ef = l;nf(^) = 100; E c 2 = 2;n c 2 (E%) = 
200. The first two simulations with both a — 0.005 but 
different random numbers in H mt show fast relaxation. 
The smaller the entries in the interaction-Hamiltonian, 
the slower equilibrium will be reached which is shown in 
the third simulation using a — 0.001. The initial states 
for these simulations are \\P(t = 0)) = ® |?i) c with 
Hq < n < Hq + n\ . 
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time 



Fig. 4. '2-level-gas' under canonical constraints: Simulations 
for three different interactions and the same initial state in g. 



The equilibrium value for this system may be calcu- 
lated using Boltzmann's a priori postulate applied to the 
energy eigenstates of the total quantum system: 



p 9 (E 9 



Poo 



: p 9 {E 9 = 0) 



n c (U - E 9 )n 9 (E 9 ) ^ 
n tot (U) : 



n c (2-0)n 9 (0) 200-1 



(6) 



_ 2 

f (2) ~ 200 + 100 ~ 3 
. (7) 

which is in surprisingly good accord with our simulations; 
the latter also show the important fact of independence 
of the actual microscopic details of the interaction, as re- 
quired for analogies to thermodynamic behaviour. This 
fact does not depend on any assumptions about the struc- 
ture of degeneracy in the system but is a general aspect 
of composite quantum systems. We have checked and con- 
firmed the independence of the equilibrium value from the 
starting level n in the container. Furthermore, neither the 
equilibrium value nor its variance depend sensitively on 
a, the parameter for the strength of the coupling, as long 
as a <C 1. We have seen deviations from thermodynamic 
behaviour for some large a, but - without supporting a 
clear-cut message - this parameter-space has not been in- 
cluded in our present investigation. 

However, the full independence of the initial condi- 
tions, i.e. of the initial energy distribution in the subsys- 
tem g can only be reached for 'thermodynamic schemes' 
of degeneracy (defined above), as can easily be shown. In 
Fig- 111 three different starting probabilities 0,0.5 and 0.9 
in state |0) s lead to the same equilibrium value of |. 

To emphasize the distinction between our model and 
common 'system-bath' models or open system models 6 , 
we calculate the occupation of the container levels in these 
simulations. Whereas the latter models take for granted 
a thermal (Boltzmann) distribution of the energy in the 
bath, we obtain (i) different distributions for different ini- 
tial conditions (but same distributions in g) and (ii) dis- 
tributions that are far from any thermal distribution (e.g. 
for the case of po(|0) 9 ) = 0.5 the equilibrium occupation 
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Fig. 5. '2-level-gas' as of fig. El Simulations for three different Fig. 7. Variances of the first diagonal element of the reduced 

density matrix of the gas over system size for 13 simulations 
plus a least square fit (line). 
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Fig. 8. Energy spectrum used for the two simulations repro- 
ducing the Boltzmann-distribution. The two dots symbolize 



Fig. 6. '2-level-gas' as of fig. 03 The von-Neumann entropy for the sharp initial energies in both the container and the gas. 
the simulations from fig.Q3 



of the degenerate container energy state with E\ — 1 is 
larger than the occupation of the lower state with Eq = 
by a factor of 3). 

For a closer connection with thermodynamics, the von- 
Neumann entropy of the subsystem g can be calculated 
according to eq. (0). This has been done for the three 
simulations from fig. [S] and is shown in fig. El As expected 
for thermodynamic systems, the entropy increases up to 
an equilibrium value that does not depend on the initial 
conditions and fluctuates around this value (cf. Jl]). The 
equilibrium value from the simulations is in accord with 
the maximal possible value for a 2 x 2 matrix p 9 with 
the equilibrium values of Pq and p\ 1 on the diagonal and 
zeros in the offdiagonal (eq. JSJl): S 9 eq th w 0.6365^- 

All the local measures considered so far exhibit fluc- 
tuations around their equilibrium values after relaxation. 
Fig. [7] shows the variances of Pq of 13 simulations that 
all produced the same equilibrium value of approximately 
|. The variances are plotted over the level of degeneracy 
n c (Ei) of the second (starting) energy level of the con- 
tainer as a measure for system size. The line in fig. is 
the least square fit with [Ap 9 m ) 2 ss ^{E^) ; i- e - the variance 
decreases with increasing system size. This is the result 



from standard thermodynamics |13j and also follows ana- 
lytically for the quantum case of bi-partite systems |15j . 

To further clarify the connection between quantum 
mechanics of small embedded systems and thermodynamic 
behaviour, the number of levels in the gas-subsystem has 
finally been increased from two to five. We restrict our- 
selves from now on to sharp initial energies, not only in 
c, but also in g, which reduces the number of energies 
that have to be considered for the dynamics to five in 
the container (Fig. (SJ. Again, it does not matter, how the 
spectrum looks like in between those levels: The conser- 
vation of energy distribution of the Schrodinger equation 
requires only those levels shown here. 

Two simulations have been carried out based on the 
two different degeneracy schemes shown in tabled 

The degrees of degeneracy obey our 'thermodynamic 
condition' as defined above: for simulation 1, n c (E c ) = 
6 • 2 E ", for simulation 2, n c (E c ) « 17 • (1.7) E °. For both 
simulations averages of the diagonal elements of the re- 
duced density matrices are calculated after relaxation and 
are shown in fig. 

The energy distributions in the gas follow a Boltzmann- 
distribution, the expected distribution for a thermody- 
namic system under canonical coupling. The least square 
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Table 1. The two degeneracy schemes used to reproduce the 
Boltzmann-distribution. 

E c 10 12 3 4 



n c (E c ) for sim.l 
n c (E c ) for sim.2 



6 12 24 48 96 
17 29 49 84 142 
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Fig. 9. '5-level-gas' under canonical constraints: Averages of 
the probabilities to find the subsystem g in the corresponding 
energy level for the two simulations together with the least 
square fits. 



fits give p3(E g ) m 0.53 • e -°- 73B9 « 0.53 • (2.1)- £9 for sim- 
ulation 1 and p9(E g ) w 0.44 • e -°- 52E9 w 0.44 ■ (1.68)- E " 
for simulation 2. 

Through the exponents of those two fits, two different 
temperatures can be introduced to describe the two differ- 
ent systems. The value of those temperatures in our model 
is determined completely by the degeneracy scheme of the 
larger subsystem c. (Under more general conditions the 
temperature will also depend on the initial energy distri- 
bution.) Again, we did not 'produce' this result by enforc- 
ing a thermal energy distribution of adequate temperature 
into the container. Only the smaller of the two subsystems 
evolves into a thermal equilibrium state, the larger one 
does not. 



5 Summary 

We have studied small quantum systems, in which a par- 
tition into two differently sized subsystems typically leads 
to a quasi-irreversible behaviour in the smaller subsystem 
(the 'gas') that can be called thermodynamic in many 
respects, even though the whole system is described by 
a pure Schrodinger-dynamics. For our specifically chosen 
initial states the 'container' does not approach a thermal 
equilibrium state. As opposed to master equation treat- 
ments that assume no entanglement at all between the 
considered system and its environment, the only reason 
for one subsystem to approach equilibrium is the entan- 
glement between the two subsystems. 

All the results obtained here are in good accordance 
with recent attempts to base thermodynamics on the mi- 



croscopic principles of quantum mechanics |Zj |H] and 
are to be seen as a generaliziation of real physical models 
described before JUj [TT] . 

Our paper may also shed new light on the recent de- 
bate about possibilities to violate thermodynamic expec- 
tations, in particular the second law of thermodynam- 
ics jT^] 53- The conditions found in our model for the 
smaller subsystem to behave thermodynamically can, if 
broken, also be used as a gateway to non-thermodynamic 
systems. E.g. a non-exponential dependence of the degen- 
eracy of the energy would lead to energy distributions 
within the 'gas' that were no longer the Boltzmann-type. 
Furthermore, for small systems the exponential decay (in 
terms of pertinent occupation probabilities) may change 
towards what would look more like a Gaussian. Further 
investigations are in progress. 

Financial support by the Deutsche Forschungsgemein- 
schaft is gratefully acknowledged. 
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